Fundamental Limits of "Ankylography" due to Dimensional Deficiency 
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Single-shot diffractive imaging of truly 3D structures suffers from a dimensional deficiency and 
does not scale. The applicability of "ankylography" is limited to objects that are small-sized in at 
least one dimension or that are essentially 2D otherwise. 



Raines et al. [l| proposes a method, dubbed "anky- 
lography" , for three-dimensional (3D) structure determi- 
nation using single-shot diffractive imaging (SSDI). But 
the conclusion without limitation that the 3D structure 
of an object is "in principle encoded into a 2D (two- 
dimensional) diffraction pattern on the Ewald sphere" 
and may be inverted by SSDI is inadequately substan- 
tiated and conceptually misleading. Here we point out 
that SSDI in general suffers from a dimensional deficiency 
in using a 2D observation with a quadratically growing 
number of degrees of freedom (NDF) to recover a 3D 
structure with a cubically growing number of unknowns, 
when the size of a genuine 3D object increases. Practi- 
cally obtainable signal-to-noise ratio (SNR) and measure- 
ment accuracy limit the applicability of "ankylography" 
to objects that are small-sized in at least one dimension 
or that are essentially 2D otherwise. 

The rate of reliable information transfer via a spatial or 
temporal channel is fundamentally limited by the chan- 
nel capacity that is determined by the NDF available 
therein and obtainable SNR Q ■ The resolving powers of 
telescopes or radio antennas and microscopes including 
SSDI are all limited in much the same manner. A steep 
(exponential) price in signal power has to be paid to ob- 
tain data rates or resolutions significantly beyond that 
are supported by the available NDF 0, y]. SSDI uses a 
spatial channel characterized by a bounded linear oper- 
ator T : L 2 ([0, 1] 3 ) — > L 2 {S) that projects a real space 
amplitude with support [0, 1] 3 onto the Ewald sphere 
S = {(/, 5, h) : f 2 + g 2 + h 2 — 1}, where I is measured in 
units of the wavelength, (/, g, h) represents a normalized 
spatial frequency. Well known is the existence of a pair 
of orthonormal bases {wi}i>i and {vi}i>\, called normal 
modes, and the associated modal gains {Ai(T)}j>i, be- 
ing all positive and arranged in a nondecreasing order, 
such that Tui — \f\ivi, \/i > 1 0, 0]- For any modal 
cutoff threshold e £ (0,1) as determined by practically 
obtainable SNR and accuracy in signal measurement, 
N(T, e) = max{ i : A»(T) > e} is the number of usable 
normal modes. We shall prove that when I is large, 

N(T, e) < 8l 2 + 0(| log e| Hog I + log 2 e log 2 I), (1) 

which grows far too slow in comparison with the number 
of unknowns 0(l 3 ) in the structure of a general 3D ob- 
ject. The dearth of NDF would persist even if 1/e grew 
exponentially as I increased, so long as the exponent grew 



no faster than I 3 / 2 , more specifically, | loge| = 0(l 15 ~ s ), 
for any fixed small <5 > 0. 

By considering the distribution of received energy be- 
tween two hemispheres, it is easy to see that N(T, e) < 
N(T+,e/2) + N(J-,e/2) < 2JV(T D ,e/2), where T± : 
L 2 ([0,Z] 3 ) -> L 2 (S ± ) and T D : L 2 ([0,l} 3 ) L 2 (S D ) 
are restrictions and extension of T to the corresponding 
codomains S± — {(/, g, h) : f 2 +g 2 + h 2 = 1, h ^ 0} and 

Sn = {(/,<?) ■■ l/l < 1, Iffl < 1} = [-1, 1] 2 - T P = T z F y F x 
is a product of three linear operators, with 

[F y F xP }(f,g,z) = I J J p(x,y,z)e- i2 ^ x +^d X dy 

being a bounded linear operator of norm 1 , and 

[T z a](f,g) = ±J^(f,g,z)e- i2 ^-f 2 -^ 1/2 -^dz 

being also a bounded linear operator of unit L 2 operator 
norm 5]. Well known results for operators of "time and 
frequency limiting" Q state that 

N(F x ,e) = N(F yi e) = 2l + 0(\ log e| log/), (2) 

and the x-y separability implies that 

N(F y F x ,e) > N(F yi e^ 2 )N(F Xl e^ 2 ), (3) 
N(F y F x ,e) < N(F y ,e)N(F x ,e), (4) 

consequently, 

N(F y F x ,e) = [21 + 0(| log e| log I)} 2 
= 4/ 2 + O(|loge|nog/ + log 2 elog 2 0. (5) 

However, the T z term makes the singular value distribu- 
tion of To more complicated. Fortunately, an operator 
inequality comes to the rescue Q- It follows from 

t*t d < ||T Z || 2 f;f;f,,f, = f;f;f„f, (6) 

that Aj(Trj) < Xi(F y F x ), Vi > 1. It is now obvious that 

N(T,e) < 2N(T n ,e/2)<2N(F y F x ,e/2) 
= 8/ 2 + 0(|loge|nog/ + log 2 elog 2 Z), (7) 

which constitutes a rigorous proof of equation ([1} . 

It is worth noting the fundamental nature of the limi- 
tation, that single-shot diffraction does not convey suffi- 
cient information to invert the 3D structure of an object, 
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even the amplitude (instead of intensity) of the diffracted 
field is sampled continuously and measured directly with 
no phase ambiguity. With practically obtainable SNR 
and measurement accuracy that determine a threshold e 
of modal cutoff, any signal in the linear space spanned 
by the normal modes of orders higher than N(T, e) is 
essentially lost in transmission or attenuated beyond de- 
tection. Oversampling and inversion algorithms are irrel- 
evant in this context. Indeed, the finiteness of an object 
ensures that the entire diffraction field is uniquely deter- 
mined by a finite number of Nyquist-sampled amplitude 
values. With a limited NDF, the only way to transmit 
and retrieve more information is to increase the number 
of significant figures in the measured amplitudes, which 
however quickly becomes prohibitively expensive. 

In summary, SSDI of truly 3D structures does not 
scale. The applicability of "ankylography" is limited to 
objects that are small-sized with respect to the wave- 
length in at least one dimension or have structures being 
essentially 2D in complexity. That may be the case in 
Raines et al.'s computer tests and preliminary experi- 
ment [lj. Raines et al. [l[ also put much emphasis on 
certain "physical constraints" , many of which as being 
described and used are actually steps of numerical pro- 
cedures instead of mathematical constraints of model for- 
mulation, while the bona fide physical constraints of con- 
tinuity and boundedness in the support and uniformity 



outside are automatically satisfied by the normal modes 
in the present formulation. The nonnegativity of the ob- 
ject field fixes only a single degree of freedom, i.e., a 
global level shift. The incorporation of more physical 
constraints, to the extreme of having the majority of the 
0(l 3 ) unknowns fixed, could arguably alleviate the prob- 
lem of dimensional deficiency, however that diminishes 
the generality and appeal of "ankylography" . 
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